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' Abstract 



We study spontaneous scalarization in Scalar- Tensor boson stars. We find that scalariza- 
' tion does not occur in stars whose bosons have no self-interaction. We introduce a quartic 

self-interaction term into the boson Lagrangian and show that when this term is large, 
scalarization does occur. Strong self-interaction leads to a large value of the compactness (or 
\ sensitivity) of the boson star, a necessary condition for scalarization to occur, and we derive 

an analytical expression for computing the sensitivity of a boson star in Brans-Dicke theory 
from its mass and particle number. Next we comment on how one can use the sensitivity 
of a star in any Scalar- Tensor theory to determine how its mass changes when it undergoes 
gravitational evolution. Finally, in the Appendix, we derive the most general form of the 
^ ■ boson wavefunction that minimises the energy of the star when the bosons carry a U{1) 

charge. 
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^ ! 1 Introduction 



Scalar Tensor (ST) theories of gravity are the most natural generalisations of General Relativity 
(GR) and describe gravity as being mediated by both a metric field gab and a scalar field In 
the simplest ST theories, Newton's constant G is replaced by a field at the level of the 

action, while $ obeys a wave equation and is coupled to the curvature via a function u;($). The 
strength of this coupling decreases as uj increases and the special case of constant to gives the 
Brans Dicke (BD) theory. One recovers GR by taking the limit u; — > oo, <I> — > ^gRi where ^gr 
is constant. By choosing appropriate units, one can set ^gr = 1. In the original (Jordan frame) 
formulation of these theories, the ordinary matter is universally coupled to the curvature, so 
these theories are metric theories. As is well known, a conformal transformation may be used 
to express the theory in a new set of variables 5*^ = ^Qab^ ^ = fi^)^ where 5*^ is the Einstein 
frame metric and (/? is a new scalar field that plays the role of a matter source and also couples 
the matter to g*^. The Einstein frame representation is often more convenient mathematically 
than the Jordan frame (for example, the Cauchy problem is well posed in the Einstein frame 
||^]) although the Jordan frame is often considered to be the physical one. 

At present, the primary motive for studying ST gravity comes from supergravity and super- 
string theories whose low energy limits always include one or more scalar fields (dilaton fields) 
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that play a similar role to the ^ field of ST gravity. However, the coupling between the non- 
dilaton terms and the curvature in string actions is non-universal so that the correspondence 
between the low energy limit of supersymmetric theories and ST gravity is only approximate. 
Nevertheless, this approximation is fairly accurate since, to date, highly sensitive experiments 
have detected no violation of the principle of universal coupling Q . 

In general, the results of Solar System experiments place the strongest constraints on the 
viability of most theories of gravity. For metric theories, these experiments are interpreted 
within the framework of the parameterised post-Newtonian (PPN) formalism which describes 
the gravitational field of the Solar System as an expansion about flat spacetime in powers of the 
Newtonian gravitational potential U near to the surface of the Sun. The coefficients appearing 
in the expansion up to the order terms, the so called first post Newtonian (IPN) order, 
are the 10 PPN parameters. In ST gravity only two of these are non-zero: 7 and /?, while 
is everywhere close to some cosmologically determined boundary value ^b- In terms of the 
coupling function to the two PPN parameters of ST theory are given by 

7 = ^, P = l + {2u.B + 3r\2u.B + ^r'^, (1) 
2 + ujb d<$ 

where ujb = uj{^b) and doj/d^ is to be evaluated at ^ = ^b- The observational constraints on 
7 and (3 are 

I7 - 1| < 0.0003, 1/3 - 1| < 0.002. (2) 

For BD theory, the first of these inequalities implies that uj > 3300 and all of the predictions of 
BD theory differ from those of GR to within a relative error of ~ 1/uj. For other ST theories, 
the observational data only places limits on the behaviour of in the slow motion, weak 

field limit. 

A framework for analysing the motion of compact relativistic bodies in metric theories of 
gravity, the modified EIH formalism, has been developed (see Eardley |^] and Will Q for details). 
The formalism treats each body as a point mass moving in the gravitational field of its companion 
bodies. Each body has an inertial mass which governs its response to the ambient gravitational 
field via quasi-Newtonian equations of motion. The role of non-metric gravitational fields (such 
as the ^> field in ST gravity) is to modify the inertial mass. The modified EIH formalism forms 
the basis of the post Keplerian formalism, which gives a parameterised description of the motion 
of binary pulsar systems and is used to interpret observational data gathered from these systems. 
These data must be used to constrain or rule out theories of gravity on a case by case basis (in 
contrast with solar system observations which can place limits on entire classes of theories by 
limiting the values of the PPN parameters). 

Since the coupling function lv is completely arbitrary, there are an infinite number of viable 
ST theories. This leaves open the possibility that gravity is described by a ST theory whose 
predictions are arbitrarily close to those of GR at the current epoch and in the weak field 
limit, but differ considerably in strong field situations or at earlier cosmological epochs. Thus 
astrophysical objects with strong internal gravitational fields may still exhibit behaviour that 
differs greatly from that predicted by GR. 

One strong field effect is spontaneous scalarization, recently discovered in neutron stars by 
Damour & Esposito-Farese ||5|, They found that, for ST theories in which the coupling 
function obeys the inequality 



~ ^ 2^B dw 
^ {2uj + 3)2 



< -4 (3) 
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and satisfies the constraints of eqn (|2|), the Einstein frame scalar field inside a neutron star 
rapidly becomes inhomogeneous once the star's mass increases above some critical value. For 
a star whose mass is below this value, is nearly constant throughout the star (a state which 
minimises the star's energy), while for higher mass stars, the energy is minimised when (p has 
a large spatial variation. These effects become more pronounced in the limit $b — > ^gr^ 
lob ^ oo. The coupling function they chose was of the form 2a; + 3 ~ l/(log ^) which diverges 
as <I> ^ ^GR = 1 a-iid gives weak field neutron stars that are indistinguishable from their 
GR counterparts. A similar study was performed by Salgado, Sudarsky & Nucamendi who 
considered neutron stars in a theory with a non-linearly coupled scalar field x which is equivalent 
to standard ST theory with uj ~ — 1), which also diverges as <l> — > ^gr- Independently of 

Damour & Esposito-Farese, a study of strong field effects in ST neutron stars was carried out 
by Zaglauer However, the formalism developed in involved an artificial definition of the 
sensitivity (defined below) of a neutron star, leading to the apparent occurrence of strong ST 
effects for choices of uj and that the study in |^ showed should not occur. 

It turns out that when spontaneous scalarization effects are present in neutron stars, the 
modified EIH formalism breaks down. An alternative, more sophisticated formalism, adapted 
exclusively to tensor-multi-scalar theories of gravity, has recently been devised by Damour & 
Esposito-Farese Q. Their formalism, which is based on an expansion of an effective coupling 
strength a between the ip field and the non-scalar matter about its weak field value, is well 
suited to the study of the scalarization phenomenon. 

The purpose of the present paper is twofold. First we discuss some aspects of the modified 
EIH formalism as it applies to ST gravity after introducing appropriate definitions of the inertial 
mass and sensitivity of a star. We then numerically compute the sensitivity of boson stars for 
several choices of to. Secondly, we demonstrate the existence of spontaneous scalarization in 
boson stars and compare the behaviour of the sensitivity and the coupling parameter a with 
their behaviour in neutron stars. 



2 Field Equations and Scalar Field Boundary Values 



The Jordan frame action for scalar-tensor gravity universally coupled to one or more matter 
fields fi is given by 



d —g 



(4) 



where gah is the Jordan frame metric, R is the Ricci scalar curvature formed from g^h and we 
use units in which G = c = 1. The generalised Einstein equations in this frame are found by 
varying / with respect to g""^ and are given by 

Gah = ^(VaVb$ - gab9"'VcVd<^) + ^{da^Ob<^ " lgab9"'dM<^) + ^T^b (5) 



where 



Tab '■- 



2 5{yf^Lm) 



(6) 



is the matter encrgy-inoniGntiini tensor, Gab the Einstein tensor formed, from g^fj and Vq, is 
the covariant derivative compatible with gab- The scalar and matter fields obey the equations 



1 



2u; + 3 



STTT-^g'^'dM'^ 



VaT"^ = 



(7) 
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where T := g°'^Tab- 

To rewrite the field equations in the Einstein frame, we make the field redefinitions 



9:, = A-^g^,^ gf = A^g^\ dv = ^^^^d^, (8) 

where 

A^{^) ■= (9) 

is found by integrating the third of eqns ^ and solving for <I>. In the Einstein frame, the 
generalised Einstein equations are 

Gl, = 2daipd,^ - glbdfdc^dd^ + SttT:, (10) 

where 

r;, := A^T^,, (11) 

while the matter and scalar field equations are 

afvyi^ = -^^T*, VlTf = -^^5.^ Tf, (12) 

where T* := g"'^T*^ and V* is the covariant derivative compatible with g*^^. 

We analyse the behaviour of boson stars for three particular coupling functions to. The first 
is the BD coupling, where to is constant. The second is the exponential coupling law chosen in 
1^], which in the Einstein frame is 

A{ip) = exp{-kip^), (13) 
where A; is a constant. This is equivalent to the Jordan frame coupling 

2u; + 3 = — ^ . (14) 

2k log ^> ^ ^ 

Our third choice is equivalent to that made by Salgado et al, who consider a theory in which 
the Lagrangian is given by 



167r 



(1 + Cx')R - Ig'^'daxdbX 



+ Lm, (15) 



where ^ is a constant and x ^ non-linearly coupled scalar field which is constrained to have the 
value X ^ 0- Making the field redefinition 1 + = ^i one finds that the theory derived from 
eqn (jl^) is the standard ST theory outlined above with the coupling function 

W = -TT T, (16) 

8^(^-1) ^ ^ 

where $ > 1. It is more natural to choose ^ instead of x the physical scalar gravitational 
field since, when studying the motion of test bodies about a gravitating central source, it is ^ 
that appears explicitly in the equations of motion for the test bodies. 

We consider asymptotically flat solutions in which $ tends towards some limiting value $oo 
at spacelike infinity. We assume that a star embedded in a cosmological background can be 
approximated by an asymptotically flat solution and we identify $oo with $b at the present 
epoch. Then equation (^) places restrictions on the maximum allowed value of <I>oo for any 
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solution at the current epoch. Using eqn (||), one finds that for the couphng given by eqn (|T^), 
the maximum allowed boundary value is given by the lower of the two values 

(f^max _ 1 _|_ 0-0003 ^max _ -i _|_ O-OO^ /-, ^\ 



Similarly, for the coupling function given by eqn (|16D, the limiting value of $00 is given by the 
lower of the two values 

^max _ 0m03 ^1+0^ (18) 

3 Mass and Scalar Charge 

To study the mass of a boson star, or any other stellar object, we may consider it either as a 
central source of gravitational field and explore the surrounding spacetime with test particles 
of differing compositions, or consider the star to be an extended, self-gravitating test body and 
examine its motion in the gravitational field of a second massive body. Either situation allows 
one to define several related quantities based upon the mass of the star and the behaviour of <^ 
in the far field region. The former situation was assumed by Salgado et al. Q , who examined the 
mass and scalar charge (defined below) of neutron stars considered as central sources. The latter 
situation is described by the modified EIH formalism and the formalism devised by Damour and 
Esposito-Farese and is considered by them in js], We shall briefly summarise deflnitions of 
mass and scalar charge for both situations, as applied to any compact object (such as a boson 
star or neutron star), and outline part of the modified EIH formalism. We specialise to systems 
consisting of spherically symmetric bodies each of whose line elements, when described in the 
Jordan frame, are given by 

ds'^ = -B(r)dt^ + ^——dr"^ + r^dS^ (19) 

^ ' l-2m{r)/r ^ ' 

where dS'^ is the line element of the unit 2-sphere and m(r) is the usual generalised Schwarzschild 
mass. Each body is modelled as a static, asymptotically fiat object. 

Consider first the motion of test particles about a central star. We define the Jordan frame 
scalar charge of the star by 



Qs := lim 



2^$ 

dr 



(20) 



Taking the ADM mass M^dm '■= lirar^oo'm{f) together with the scalar charge, we define the 
tensor mass 

Mt := Madm - ^ (21) 

where $00 '■= lim^^oo ^ is the asymptotic value of the gravitational scalar field. Mt measures 
the total energy of the star and is also the active gravitational mass measured by an orbiting 
test particle for which $ is locally homogeneous (such as a small black hole). An orbiting test 
particle with negligible self-energy moves on geodesies of gab and measures the Keplerian mass 

Mk := Madm - (22) 

This is different from Mt since the Strong Equivalence Principle (SEP, defined in Q) is violated 
in ST gravity. 
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From the expressions for A'lx and Mt we have 



Qs = 2<^^{Mt-Mk). 



(23) 



A locaUy freely fahing Cavendish experiment far from the star measures a gravitational coupling 
strength 

.1 2uJoo + 4 



G = $ 



°° 2woo + 3 



where cooo = a;($oo)- Given a particular form of u, one can then calculate $c 
combining the above two equations, we find that Qs is measurable in principle. 
Salgado et al. define and use the scalar charge 



lim 



dr 



(24) 
Hence by 

(25) 



associated with the field x that appears in the Lagrangian ([T5|). This quantity is related to the 
charge Qs associated with <I> by 



Qs = 2ixQs = 2Qs^/Wc 



1 . 



(26) 



From the above definitions it appears more natural to use Qs as the physical scalar charge of 
the star instead of Qs- 

Both Mt and Qs may be related to corresponding Einstein frame quantities as follows. In 



the Einstein frame, the Schwarzschild radial coordinate r* is given by r* 
defining the Einstein frame scalar charge Q*^ by 



Q\ 



lim 



dr* 



and using the third of eqns (P) gives 



Q*s 



Qs- 



One can also show that Mt is related to the Einstein frame ADM mass M'\j~ij^ by 



M 



ADM 



'^ooMt- 



Then, 



(27) 



(28) 



(29) 



Consider now the motion of a spherical star, with some conserved particle number N ^ in the 
gravitational field of a companion body. We assume that the separation between the two bodies 
is much greater that either of their radii, and we denote the inter body value of the scalar field 
by $oo- The scalar field in the interior of the star must be matched smoothly to $oo and 
its value will affect the structure of the star. In particular, its inertial mass M/ and hence its 
motion in the gravitational field of the companion body will be sensitive to changes in the value 
of <I>oo- Expanding about as = <I>oo + 5<5, we write 



Mi = Mj^ + ?^<5$ + 0(5$2) 



where Mj^ is the value of Mj when ^t '- 

Mi = Mi^ 



9$ 



^^oo, or 



1 + s 



(5$ 



+ 0(5$^ 



(30) 



(31) 
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where s is the first sensitivity defined by 



_ djlog Mj) 

d{log<^) ■ ^'^ ' 

In eqns ( pO| ) and (p2|), the derivatives are evaluated at $ = $00 for fixed A'^. Higher order 
sensitivities may be defined by taking successively higher order derivatives of Mj, to obtain the 
expansion for Mj to higher powers of 5^. s measures the response of a star's inertial mass to 
changes in the value of $00 and it also gives a useful measure of the compactness of the star 
Here we define M/ to be 

Mj := $ooMt. (33) 

This differs from the definition given by Will |^] by a factor of {2uiao + 3)/(2a;oo + 4) and is 
chosen so that our definition corresponds with the Einstein frame quantity defined by Damour 
& Esposito-Farese in Q. The difference in definitions of Mj does not affect the BD analysis 
given in |Q] since, for constant uj, this factor cancels in the definition of s given there. For more 
general ST theories, our definition if Mj leads to a self consistent set of equations of motion 
that, to IPK order, are identical to the BD equations derived in Q (except for the different 
definition of M/) and obviates the need to introduce the modification of s given in 

These definitions form the basis of the modified EIH formalism, examined in detail for more 
general theories of gravity in Q , and s and its derivatives play a fundamental role in determining 
the motion of a body in the gravitational field of other bodies. Formally, to derive the equations 
of motion governing such a system, one first determines the inter body gravitational fields from 
the field equations where the Mj appear as delta function terms in Tab- One then substitutes 
the solution into a quasi-Newtonian particle Lagrangian, taking into account changes in the Mj 
due to the relative velocities of the bodies, to find the equations of motion for each body. These 
equations explicitly involve the sensitivities of the bodies, which determine how the motion 
deviates from the Newtonian prediction, and are written as an expansion in powers of f /c about 
the Newtonian solution, where u is a characteristic velocity of the system. (One cannot expand 
in powers of the Newtonian potential of the bodies since for strongly bound objects, U is of 
the order of 1). The expansion is normally taken to order v'^/c^, the first post Keplerian (IPK) 
order, and in the limit that U for each body becomes small, the equations of motion reduce to 
the IPN equations. The post Keplerian equations then are a parameterised description of the 
modified EIH equations of motion, adapted so that the system is described directly in terms of 
observable parameters. 

For a system composed of two bodies it turns out that only the first sensitivity s appears 
in the equations of motion to IPK order. For example, as shown in detail in [^, given a 
star of inertial mass Mf and sensitivity sa with a companion of mass M^ and sensitivity sb 
at a distance Rab, the quasi-Newtonian gravitational force F between the two bodies has a 
magnitude F = GAsMfM^ /R\g where the effective gravitational coupling strength Gab is 



^ 2a;oo + 4 

Gab 



$oo(2a;oo +3) 



1 , 

1 -^(■5A + SB- 2SASB) 



(34) 



to IPK order. (Note that the above equation is a simplification that ignores effects due to the 
relative velocity of the two bodies.) If body _B is a non-self-gravitating test particle then sb = 
and body A has a Kepler measured mass 



^A ^ FRab ^ (2^00 + 4)M/ / _ SA \ 
^ Mf $oo(2a;oo+3) V oj^ + ^J 



(35) 
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In ST gravity, the $ field outside an asymptotically flat static black hole is constant and one 
can show that the sensitivity s bh of a black hole is s bh = \ ■ Assuming that a small black hole 
in orbit about body A also has this property, then the orbital mass in this case is given by 

(2u;oo + 4)M/ / _ 1 \ ^ Ml_ 



^ . - , (36) 

$oo(2woo + 3) V 2a;oo+47 $c 



which consistent with our relation (|3^) between the inertial and tensor masses. 

The formalism devised by Damour & Esposito-Farese Q for tensor-multi-scalar theories is 
similar to the modified EIH formalism except that it describes the motion of the system in terms 
of Einstein frame variables. In place of s, the fundamental quantity used is a coupling parameter 
a which, for a star of mass M^j-^j^j and Einstein frame scalar charge Q^, is defined by 

and measures the strength of the coupling between ip and the non-scalar matter. In the weak 
field limit, a — > — l/\/2a;oo + 3. The equations of motion at the IPK level then depend upon 
the values of a and da/ dip ^ evaluated at = (^oo- 

Damour & Esposito-Farese (see Appendix A of Q) have shown that, by varying the Einstein 
frame Hamiltonian with respect to changes in the asymptotic value of the scalar field (p, one can 
derive the following relation: 

nt. 1 BM*^,. 

(38) 



QS _ 1 S^ADM 



where the derivative is evaluated at 93 = ^p^ for fixed particle number. A similar analysis, 
performed on the Jordan frame Hamiltonian, gives 

Qs -: 



Mi 2looo + 3 



(l-2s). (39) 



Using eqn ( |39| ) one can calculate s directly from the scalar charge and inertial mass. Combining 
eqns (|38|) and ( |39| ) and using eqns (p9|) and (|33| ) , we find that the parameters a and s are related 
by 

a = -==i==(l-2s), (40) 
Viwoo + 3 

to IPK order. 

Formally ST gravity reduces to GR in the limit u) ^ 00, $ — > ^gr = 1 (in units with G = 1 
and provided g"'^Tab is non-zero). In this limit, s generally approaches a non-zero value and 
measures the response of M/ to a (global) change in the value of Newton's constant G. 

The modified EIH formalism is adequate for describing the motion of compact bodies when 
effects such as spontaneous scalarization are absent. However, as shown in ||5|, ^, when scalariza- 
tion occurs in a neutron star, both a and Q*g remain finite and non-zero as a; 00 — > 00. Equation 
( [ioD then implies that s must diverge and so the modified EIH formalism breaks down. We 
shall discuss this further in Section |5|, where we also show how s behaves for various boson star 
models. 



4 Boson Stars 

Boson stars are gravitationally bound configurations of zero temperature bosons and many of 
their properties in ST gravity have already been studied [0, 11, O, 0|. They share many 
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features in common with simple neutron star and white dwarf solutions. Sequences of boson 
star solutions may be parameterised by the central boson field density, and when plotted against 
this parameter, the mass and conserved charge curves show a maximum at which instability first 
occurs. 

We take as a matter source a complex self-interacting boson field ^ with Lagrangian 



Lr 



'2 



(41) 



where A measures the strength of the boson self-interaction and we have chosen units in which 
h/ ^ = 1, where /i is the boson mass. The Jordan frame energy momentum tensor for this field 
is 



Tab = 9a^a,^ + da^db^ - IgabQ^^d,^ + dc^Od^) - 5afe*^ " ^T^QabK"^^? (42) 

We consider static, spherical solutions with the line element (|l^) and we only consider 
minimum energy solutions, which implies that the boson field wave function may be written as 



P{r) 



exp(ir2t) 



(43) 



(see Appendix) where P{r) is a real dimensionless function and the constant is real. For the 
metric given implicitly in the line element (|l9|), the independent components of the field and 
matter equations (^) and (0) are 



m 



2u; + 3 



(r — 2m) 



$'2 doj P'2 A rP2 



l + {2u + 5) 



B 



+{2oj - 1 



+ (r — 2m) 



r^'B' ^>'2rw 

+ 



AB^ 



4$2 



(44) 



B' 



1 



2$ + r$' 



^''^ujBr 



rm 



m 



r(r — 2m) 



45$' + 2P'^Br 

B'A 2 
~2B~ 2UJ + 2, 



2m 



2P'^r'^(B -n^) + P''BA + 



4Bm^ 



^>'2 duj ,n P^r 

\- P + - 

2 d$ r 



and 



P" = P' 



rm 



m 



r{r — 2m) r 



B' 
2B 



+ 



Pr 



r — 2m 



2m 
B 



— + 2AP^ 
B 



+ Ap2 



(45) 
(46) 

(47) 



where a prime denotes d/dr. For the corresponding Einstein frame equations, see |1 1[ ]. 

The existence of a global U{1) symmetry in the matter Lagrangian (41) implies the existence 
of a conserved charge N. Using the above coordinates and field variables, this may be written 
as 

r^np^ 



N 



dr 



(48) 



VB^/I - 2m/r 

and is interpreted as the total number of bosons in the star. 

must be integrated numerically. To ensure that the solutions describe 



Equations (44) to 
bound eigenstates of ^ one must impose the boundary condition 



lim P{r) 

1 — >oo 







(49) 
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and the regularity conditions 

mo = 0, P^ = 0, $^ = 0, (50) 

where the subscript 'o' denotes values at the origin. The field equations then become eigenvalue 
equations for Q, are parameterised by Pq for fixed ^>oo and automatically lead to the boundary 
conditions $ = $00 + 0{r~^), m{r) = Madm + 0{r~^) and B = B^o + 0{r~^) as r ^ cxd. For 
any solution, one can make the rescaling B^ — > 1 by rescaling 17 appropriately. 

In Brans-Dicke theory, where dco/d^ = 0, the equations with A = are invariant under the 
rescaling 

kP, ^> ^ (51) 

where k is a constant. Equation (|5^) leaves Mt, Mk and Madm invariant and rescales the 
particle number as 

N K^N. (52) 

We can use this scaling invariance to find an alternative expression of s in a A = BD boson 
star as follows. Consider a pair of solutions, cJi and cr2, with the same boundary value and 
with central boson field amplitudes Pg^'* and P^^ = Pq^"* + 5Pq, such that 5Po <C 1. Solution 
cJi has mass M^^ and particle number N^^\ while C72 has mass m!^^ and particle number A^^-^^. 
Then, to first order in 5Pq, 

= M? + ^5Po (53) 

and 

Ar(2) = + l^jPg, (54) 

where the derivatives are taken with <I>oo held fixed. We use eqn (51) to generate a new solution 
a2 with mass particle number A^*^^) and boundary scalar field $00 and we choose a value 

of K such that iV*-^) := k^N^"^^ = N^^\ From eqn (^), must then be given by 

to first order in 5Pq. This implies that 

2. ^ ^ dN 



^oo = At'^oo = $oo(^l-^^5PoJ, (56) 

while Mj^ is invariant under this rescaling, thus Mj. = Mj. ' . Then, using eqn (|53|), we have 

^(2) _ ,^(1) _ dMr^ 
dPo 



6Mt := M^'> - M^'> = ^SPo (57) 



and, from eqn (^), 

^$oo:=$oo-$oo = -|^|^5Po. (58) 



Combining these results with the definition ( |3^ ) gives 

^00 ON 



S'^oo (^^00 ^ °° 9Po ° 



A^(i) dPo 



6P 



(59) 
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^ " dN 
Taking the limit 6^00 — *■ and using eqn (^) we have 



(60) 



Ml \ dN J Mt ON ' ^ ^ 



where we have dropped the label from and Mt since this relation is true for any initial solution 

(62) 



o"i. One can show |12| that for fixed 



dN $00' 

Combining this result with eqn (|6^) gives an alternative expression for s: 

QN , , 

s = 1 - 63 

Note that eqns (^Tj) and ( |63D hold for any value of uj. 

When A is non-zero but finite, one cannot drive a similar result since one also needs to 
rescale A to keep the field equations invariant. However, Gunderson & Jensen |^ have shown 
that in BD theory, even a relatively small value of A causes the terms that are quartic in P to 
dominate the energy momentum tensor. For A > 100 the solutions are well approximated by 
taking the limit A ^ 00, and one can derive an expression for the sensitivity of a boson star in 
this limit. One first rewrites the field equations by making the field and coordinate re-definitions 
n = \/A-P, P = vj^Jli and then takes the limit A — > 00. The resulting equations are identical 
to eqns (|4^ ) to ( p7[ ) except that P and r are replaced by IT and p, A = 1, all terms including a 
factor of n' vanish and 11 is given by the algebraic expression 

= ^ - 1. (64) 

See Q for details in BD theory. The equations defining N ^ Qs and Mt are identical to the 
finite A equations, except that P and r are everywhere replaced by 11 and p. Then masses are 
measured in units of M'^J p? and in BD theory the new field equations are invariant under the 
rescaling 

$ K^^>, p Kp, (65) 



where k is constant. Under eqn (|6q), the mass and particle number rescale as 

Mt kMt, N ^ k^N. (66) 



Performing a similar analysis to the one given above for the A = case and using eqns (|6^) and 
), one can show that the sensitivity in the limit A ^ 00 is given by 



3 / _ J^dMT\ ^ 3 / _ Nn \ 
2 V Mt dN J 2 V ^-ooMr / ' 

The factor of 3/2 appears in these equations for two reasons: firstly, in the present case, delta^ 
differs from eqn (^8|) by a factor of 2 /3 (since the particle number N* scales differently from A'^* 
by a factor of k^/^) and, secondly, one must include an extra term in 5M^ to account for the 
rescaling of M^ in eqn (|6^), a term which does not appear for the A = case. 
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The first of these equalities is identical to the result derived for neutron stars ^ (in the latter 
case, N must be replaced by the neutron star's conserved baryon number). This is because both 
neutron stars and boson stars in the limit A — > oo have masses that scale as M^JiJ?, where for the 
neutron star /i is the baryon mass, in contrast, the mass of a boson star with no self-interaction 
scales as M^Jfi, so there is no factor of 3/2 in eqns (61) and (63). 

One cannot carry out the same analysis for more general ST theories since the ST versions of 
eqns (|5l| ) and (|65|) require uj{^) to be held fixed. Since $ rescales under eqn (]5T|), the functional 
form of u) must change to compensate in a way that depends upon the precise choice of lo. 
However, eqns (^) and (|6^) will still be approximately correct for solutions in which div/d^ 
remains small. 



5 Equilibrium Solutions 



We have numerically integrated the field equations (44) to ( p7[ ) for BD boson stars and for ST 
boson stars with the coupling functions given by eqns ([l^ ) and (16). A feature shared by all 
solutions sets it that both Mt and A'^ increase from zero at Pq = to reach coinciding maxima 
M^f^"""^ and AtC'""^) at some value of Pq = p^j"'"^). For Pq > p('"'^^), the solutions become 
unstable. The magnitude of the scalar charge \Qs\ also increases from zero at Pq = but its 
first maximum \Qs\^'^°'^^ does not occur at the same value of Pq as the maxima in N and Mt- 
For all solutions discussed here, we have found that s < 0.5 which implies that Qs, o: and 
are all negative. 

We consider first boson stars in BD theory. One can show that in the weak field limit, where 
the (dimensionless) parameter Pq ^ 0, the sensitivity is given by 



s = -2£ + 0{P^), 
where £ is the fractional binding energy defined by as 

$ooMr - N 



N 



(68) 



(69) 



For stable stars, if < and hence all stable BD stars have s > in the weak field limit. 
Numerical calculations then show that s increases with Pq to reach some maximum at a value 



of Pq > Pq 

-){max) 



(max) 



Values of s 



(max) • 



the sensitivity of the maximum mass stable solution at 



Pq = Pq , along with corresponding values a(max) of the coupling parameter a, are shown in 
Table |^ for several choices of uj for both A = and in the limit A 



oo. 



A ^ oo 


A = 




■5 (max) 


^(max) 




^(max) 




1 


0.211 


-0.258 


1 


0.144 


-0.318 


10 


0.189 


-0.130 


10 


0.126 


-0.156 


500 


0.182 


-0.020 


500 


0.121 


-0.024 


3300 


0.178 


-0.008 


3300 


0.119 


-0.009 



Table 1: Sensitivity s and coupling parameter a of maximum mass boson stars in BD gravity 
for A = and in the limit A ^ oo. 



As eqns (^) and (|6^) suggest, for each uj the sensitivity of a maximal mass boson star with 
A — > CO is greater than the sensitivity of the corresponding A = star by a factor of 3/2. This 



is also true for values of N less that the maximum value. For both choices of A, s 



(max) 



decreases 



12 



with increasing w towards some non-zero limit as a; ^ oo and to the level of accuracy quoted 
in the Table, S(^max) for a boson star in GR has the same value as for a to = 3300 star. For 
comparison, maximum mass neutron stars have been found to have sensitivities that vary from 



'(max) 



0.2 to s 



(max) 



0.39, depending upon the equation of state chosen 



Hence for 



boson stars with no self-interaction, the response of Mj to changes in ^ is much smaller than for 
neutron stars, while for boson stars with A ^ oo, the sensitivity is almost comparable to that 
of a neutron star. Note that, from the relation (pO|), we have a ^ as w ^ oo since s remains 
finite in this limit for all BD solutions. 

We next consider boson stars in more general ST theories. One can show that, for any ST 
theory and for any value of A, in the weak field limit s is given by 



-£:(2 + (2c^oo + 3)/?) + 0(Po')- 



(70) 



From eqns (|l|) and the limits (|2|), we have for any weak field solution s < when (3 < —0.0003. 

Figure || shows sensitivity s against particle number A'^ for several sets of boson star solutions 
with A = 0, for the coupling function ( p^ ) with the parameter choices k = 1 and k = 3. For this 
coupling function, /3 = —2k and the inequality (^) is satisfied when k > 2. For the values of k 
we are using, it turns out that the first inequality of eqns ([l^) places the tightest constraints 
on the value of <I>oo and these are (^oo - 1) < 7.5 X 10"^ for A; = 1 and ($oo - 1) < 8.3 x 10^*^ 
for k = 3. This latter value corresponds to the limit (p^o < 0.0012 which is smaller that the 
limit ipoD ^ 0.0043 quoted in since the constraint on the PPN parameter 7 has recently been 
tightened Q. For the sake of comparison, in the Figure we have included k = 3 solutions with 
the same boundary value used in |^ , although we also choose values compatible with the current 
observational constraints. 

The curves in Figure § show the coupling parameter a against A^ for the k = 3 solutions 
shown in Figure |l[ In contrast with the corresponding neutron stars curves shown in |5|, the 
maximum value of —a occurs when Pq > p^^™'^^) or, equivalently, when N < _/v(™<*^) and the 
solutions are unstable. This is because, for any set of boson star solutions, {Qsl^"^""^^ occurs at 



a parameter value Pq > Pq 



(max) 



k 


P 


$00-1 


^00 


^00 


^(max) 


^(max) 


^(max) 


^*)(max) 


1 


-2 


7.5 X 10"^ 


4.3 X 10~^ 


6759 


-0.0035 


-0.00920 


1.0 X 10' 


-4 


0.0058 


3 


-6 


1.1 X 10-^ 


4.3 X 10-'^ 


750 


-2.70 


-0.165 


5.3 X 10" 


-3 


0.10 


3 


-6 


1.1 X 10"*^ 


4.3 X 10"^ 


75114 


-4.00 


-0.0232 


7.6 X 10" 


^5 


0.015 


3 


-6 


1.1 X 10"^ 


4.3 X 10"^ 


7.51 X 10^ 


-4.02 


-0.00233 


7.6 X 10" 


-7 


0.0015 


3 


-6 


1.1 X 10-1° 


4.3 X 10~6 


7.51 X 10^ 


-4.02 


-0.000233 


7.6 X 10- 


-9 


0.00015 



Table 2: Numerical data for boson stars with A 
denotes values at the maximum mass solution. 



and 2tj + 3 = I /{2k log The label (max) 



Some of the data taken from the numerical calculations are shown in Table ^. Together 
with Figures || and ^, they show that the scalarization phenomenon does not occur in A = 
boson stars. As the study in ||5| shows, in the limit {^00 — 1) — > and for neutron stars with 
P < —4, the values of both —a and Qg for a star of given baryon number approach some finite 
limit (which is either zero for stars whose mass is below some critical value or non-zero for stars 



whose mass is above this value). From eqn (40), this implies that, for stars above the critical 
mass, s is negative and also that s diverges as ^/uJ^ in this limit. In addition, from eqn (|2^), 
Qs vanishes as l/^/uJ^. In contrast, for boson stars without self-interaction we find that for 
all solutions s reaches some finite, non-zero limit as (^>oo — 1) ^ 0. This implies that a must 
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vanish in this hmit for all solutions, as shown by eqn (|40|). As shown in Table ^, the magnitude 

of both Qg^"^^ and Q^^™""^^ decrease to zero as <l>oo — 1^0 and this is true in these solutions 
for all other values of A^. 

We next consider A = boson stars with the coupling function (|l6|). Assuming that — 1 
is small, we have /3 = —4,^ which implies that the inequality (|^ is satisfied when ^ > 1. Figures 
^ and ^ show curves of s and a against A^ for solutions with ^ = 1 and ^ = 2. For these 
choices of ^, the first of eqns (|T8| ) places the tightest limits on the values of ^oo-, which are 
($oo - 1) < 3.75 X 10-5 for ^ = 1 and ($oo - 1) < 1-87 x 10"^ for ^ = 2. In the Figures we 
have chosen boundary values compatible with these limits. The behaviour of the solutions is 
qualitatively similar to the behaviour of the solutions shown in Figures |l| and ^ and, again, the 
maximum value of —a occurs after A^ has reached its first maximum. Data taken from these 
calculations is shown in Table |3|. 





P 


^oo-l 




^(max) 








1 


-4 


3.75 X 10" 


-5 


3329 


-0.433 


-0.0229 


0.00035 


0.014 


1 


-4 


1.87 X 10" 


-5 


6679 


-0.438 


-0.0162 


0.00018 


0.010 


1 


-4 


1.91 X 10" 


-6 


65536 


-0.422 


-0.0051 


0.000018 


0.0032 


2 


-8 


1.87 X 10" 


-5 


3339 


-7.39 


-0.193 


0.0028 


0.113 


2 


-8 


9.41 X 10" 


-6 


6637 


-10.1 


-0.184 


0.0019 


0.109 


2 


-8 


2.38 X 10" 


-7 


262144 


-45.7 


-0.128 


0.00021 


0.077 



Table 3: Numerical data for boson stars with A = and w = <I>/(8^(<I> — 1)). 



Finally, we consider solutions to the field equations in the limit A ^ oo for the coupling 
function (14) and for k = 2>. As mentioned above, solutions with this limiting value of A are a 
good approximation of solutions having large but finite values of A in BD theory, and we expect 
this to be true for other ST theories. Figure ^ shows curves of a against A^ for these solutions 
and data taken from these integrations is shown in Table ^. In contrast with the A = solutions, 
these solutions show the same scalarization phenomena as the neutron stars studied in ||5|, |^. 
When the boundary value ($oo — 1) > 0, the coupling parameter a is small and negative for 
small A^ and decreases smoothly as A^ increases to its first maximum. For these solutions, Woo 



is finite and all have finite, but negative, sensitivities. In the limit (<I>c 



1) 



0, Wc 



oo 



the 



solutions divide into two classes: for those with A^ below a critical value Nc = 0.200, a vanishes 
while s decreases smoothly with increasing A^ to diverge at A^ = A'c- For stars with N > Nc, 
a is non-zero and increases rapidly with A^, while s remains divergent. The star with particle 
number Nc mass Mt = 0.189. Note that Figure ^ shows that the transition point between small 
and large values of —a becomes sharper as ($oo — 1) decreases, just as is the case for neutron 
stars. 



$oo-l 




^{max) 


0!-(max) 






1.1 X 10-4 


750 


-6.96 


-0.385 


3.51 X 10-^ 


0.068 


1.1 X 10"^ 


7.51 X 10^ 


-141 


-0.073 


7.09 X 10-6 


0.014 


1.1 X 10-1° 


7.51 X 10^ 


-1125 


-0.058 


5.66 X 10-^ 


0.011 


1.1 X 10-12 


7.51 X 10^° 


-10829 


-0.056 


5.45 X 10-^ 


0.011 



Table 4: Numerical data for boson stars in the limit A 
for k = 3. 



oo with 2uJ + 3 = l/(2A;log$) and 
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6 Discussion 



We have analysed spontaneous scalarization is ST bosons stars for several choices of the cou- 
pling function uj. We have found that scalarization does not occur when the bosons have no 
self-interaction since, in this case, the sensitivity (or compactness) of the stars is small. With 
the inclusion of a large quartic self-interaction, the stars become much more compact and spon- 
taneous scalarization occurs. The fact that this phenomenon occurs for boson stars as well as 
for simple neutron star models suggests that scalarization may be a universal characteristic of 
ST gravity. 

We have also given a brief introduction to the modified EIH formalism, in which the sensi- 
tivity plays an important role, and shown how one can calculate the sensitivity of a boson star 
(or any other compact object) in ST theory from its mass and scalar charge. As well as giving 
an indication of the compactness of a star, the sensitivity also measures the response of the 
star's inertial mass to changes in the asymptotic value of ^. Hence some of our results may also 
be applied to the study of gravitational evolution, in which the structure, and in particular the 
mass, of a star embedded in a cosmological background is forced to evolve as the value of the 
cosmological scalar field changes with cosmological time. This phenomenon was first studied for 
ST boson stars by Comer & Shinkai |11] and for BD boson stars by Torres, Schunck & Liddle 
These authors assumed that the star evolved quasi-statically and could be modelled as 
a sequence of asymptotically flat equilibrium solutions of constant particle number whose 
boundary value <^oo matches the cosmological value of $ (which in general is an increasing func- 
tion of cosmological time). These results were extended in |jT^. Gravitational evolution has 
also been shown to affect the cooling rate of white dwarfs, which provides a new method of 
constraining ST theories of gravity [p^]. 

Given a value of s for a star one immediately knows how its inertial mass will evolve, given 
the assumptions made in [0, 14, For example, in BD theory we have found that s is positive 
for all boson star solutions, and is positive for all BD white dwarf and neutron star solutions 
that we are aware of. Hence M/ will be an increasing function of cosmological time for all of 
these objects. Since M/ = $oo-^T5 the fractional binding energy equation (]69| ) shows that a BD 
boson star that is stable at the current epoch will have been stable at all earlier epochs, when 
Mi would have had a smaller value. 

For other ST theories s may be negative so that Mj will be a decreasing function of time. 
This was indeed the case for the boson stars studied by Comer & Shinkai in [ll], who modelled 
the evolution in the Einstein frame using a coupling function equivalent to eqn (13|). As we have 
seen, s < for this coupling and the mass M/ of the star will be larger at earlier epochs. As 
pointed out in |11], this implies that the fractional binding energy increases as we move further 
into the past and will eventually become positive, at which point the star becomes unstable. 
Thus, knowing the sign of s for a boson star, one can immediately see if can exist as a stable 
object at earlier cosmological epochs. These comments apply equally well to other stellar objects. 

The tensor mass (as we have defined it here) will in general decrease with time for all ST 
theories. From the definitions (|3^ ) and (^) we have 

where the derivative is evaluated at $ = $oo- For a black hole s = 0.5 and, since a black hole 
is the most compact of all stellar objects, one would expect other bodies to have s < 0.5. This 
is certainly the case for the boson star solutions discussed here. Then eqn (^) implies that Mt 
will always be a decreasing function of time in any cosmological model in which $ increases with 
cosmological time. 
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A Appendix: Minimum Energy Solutions 



A proof that eqn ( |43D leads to minimum energy solutions for static, spherically symmetric boson 
stars can be found in However, as noted by Jetzer no corresponding proof has been 



given for the case of a boson field with U{1) charge. For completeness, we generalise the proof 
in [16 1 to include charge and to allow for arbitrary static, asymptotically flat spacetimes. 



We start with the matter Lagrangian for bosons carrying charge e, which reads 

Lm = -y^' {WDb^ + Da^DhW) -m- V{^^) - (72) 

where 

Fab = dbAa - daAb (73) 

is the Faraday tensor with Aa the vector potential. 

Da := da + ieAa (74) 

is a covariant derivative operator, V is the boson self-interaction term (which we assume has 
no explicit time dependence but is otherwise left arbitrary) and the over bar denotes complex 
conjugation. We decompose the spacetime into spacelike slices S orthogonal to a timelike Killing 
vector field and write the line element as 

ds^ = -B{x')dt^ + hij{x')dx'dx^ (75) 

where i,j label spatial indices, are coordinates on E and 4" has components ^'^ = (1,0,0,0). 
One can choose a gauge in which 

Aa = {Aoix'), 0,0,0) (76) 

and one can show that, for the Lagrangian ([7^), the conserved particle number may be written 

as 

N 



f I h 

\ - [i{^dt^ - '^dt^) + 2eAo^^] d^x (77) 

where the subscript 'o' denotes the time component and h =Det(/iij). 
The boson field ^ has generalised momentum and velocity 

while the momentum and velocity associated with ^ are p and dtq. Since we are assuming that 
the spacetime is static, the Hamiltonian H for the system is given by 

H = d^x [pdtq + pdtq) -^+Hs (79) 

where I is given by eqn ^ using the matter Lagrangian (^) and Hg is a surface term that is 
independent of ^ and ^ . In an asymptotically flat spacetime, Mt = H<^^ . 

We consider a sequence of solutions to the fleld equations (|5|), with Tab derived from the 
Lagrangian (^), in which the flelds Qab, ^, ^ and ^' are held flxed while dt^ and dt^ are 
allowed to vary. Equivalently, we vary p and p while keeping the conjugate variables q and q 
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and the gravitational fields fixed. The desired form of these time derivatives is then the one that 
minimises the Hamiltonian, subject to the constraint that is conserved. Formally, we have 

6{H - QN) = (80) 

under variations 5{dt^) and 6{dt^), where is a Lagrange multiplier. Since all other fields are 
being held fixed, we only need to consider terms in H and N that involve time derivatives of 
the boson fields. These terms are 

Jt, V -d 

(81) 

= [ d^xJ^dt^dt'^ (82) 

and 

Nif,= [ d^xJ^i{^dt^ -^dt-^), (83) 

while eqn (80) is equivalent to 

5{Hm - OiVvt) = 0. (84) 
Substituting in the explicit forms of H^, and N^,, eqn (|8^ becomes 



/ ^^^V 4 [^idt'H) (5*^- + iO^) + 5{dt^) {dt^ - iO^)] = (85) 

which can only be satisfied if the equations 

dt^ - in^ = 0, dt^ + iO^ = (86) 

hold. This implies that 

^ = f{x') exp(J7t) (87) 

where / is a real function. This result holds for arbitrary V (provided V has no explicit time 
dependence) and is independent of the form of the gravitational sector of the Lagrangian. We 
note here that one can also derive this result by including the variations 5^ , 5{di'^) and their 
complex conjugates by assuming the wave equations for the boson fields are satisfied. 
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N 



Figure 1: Sensitivity s against particle number N for ST boson stars with 2a; + 3 = l/{2k log $) 
and A = 0. The curves are labelled by the value of ($oo — 1) and the solid portion of each 
represents stable solutions. The uppermost curve corresponds to the parameter choice k = 1. 
The remaining three curves show solutions with A; = 3 and, at the level of detail resolvable from 
the Figure, solutions with A; = 3 and ($oo ~ 1) < 1-1 x 10~^ are indistinguishable from the lower 
most curve shown. 
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Figure 2: Coupling parameter a against particle number N for ST boson stars with 2uj + 3 = 
l/(2felog$) and A = 0. The curves are labelled by the value of (<I>oo — 1) and all are for k = 3. 
The solid portion of each curve represents stable solutions. As (<I>oo — 1) 0, —a — for all 
solutions. 
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Figure 3: Sensitivity s against particle number N for ST boson stars with u = $/(8^(<l> — 
1)) and A = 0. Curves are labelled by the value of ($oo — 1) and the solid portion of each 

represents stable solutions. The uppermost curve corresponds to the parameter choice ^ = 1 
and is indistinguishable from ^ = 1 curves with ($(x, — 1) < 3.75 x 10~^ at the level of resolution 
in the Figure. The remaining three curves are for the parameter choice ^ = 2. 
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Figure 4: Coupling parameter a against particle number N for ST boson stars with u = 
$/(8^($ — 1)) and A = 0. The curves are labelled by the value of (<I>oo — 1) and the solid 
portion of each represent stable solutions. The lower three curves are for the parameter choice 
^ = 1 while the remaining three are for ^ = 2. 
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Figure 5: Coupling parameter a against particle number AT for boson stars with 2a; + 3 = 
l/(2/clog$) and in the limit A oo. All curves are for k = 3 and are labelled by the value of 
($oo ~ !)• The solid portion of each curve represents stable solutions. 
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